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We describe and analyze a tight-binding model of single molecule magnets (SMMs) that captures both the 
spin and spatial aspects of the SMM electronic structure. The model generalizes extended Hiickel theory to 
include the effects of spin polarization and spin-orbit coupling. For neutral and negatively charged Mni2 SMMs 
with acetate or benzoate ligands the model yields the total SMM spin, the spins of the individual Mn ions, the 
magnetic easy axis orientation, the size of the magnetic anisotropy barrier and the size of the HOMO-LUMO 
gap consistent with experiment. For neutral molecules the predicted spins and spatial locations of the HOMO are 
consistent with the results of density functional calculations. For the total spin and location of the LUMO density 
functional theory-based calculations yield varied results while the present model yield results consistent with 
experiments on negatively charged molecules. For Mni2 SMMs with thiolate- and methylsulphide-terminated 
benzoate ligands (Mni2-Ph-Th) we find the HOMO to be located on the magnetic core of the molecule, but 
(unlike for the Mni2 SMMs that have previously been studied theoretically), we predict the LUMO and near 
LUMO orbitals of Mni2-Ph-Th to be located on ligands. Therefore we predict that for these Mni2 SMMs 
resonant and off-resonant coherent transport via near-LUMO orbitals, not subject to Coulomb blockade, should 
occur. We propose that this effect can be used to identify specific experimentally realized SMM transistors 
in which the easy axis and magnetic moment are approximately parallel to the direction of the current flow. 
We also predict effective spin filtering by these SMMs to occur at low bias whether the transport is mediated 
by the HOMO that is on the magnetic core of the SMM or by near LUMO orbitals located on the nominally 
non-magnetic ligands. 

PACS numbers: 73.63.-b, 75.50.Xx, 85.75.-d, 73.23.Hk 



I. INTRODUCTION 

Single molecule magnets (SMMs) are magnetic molecules 
that contain transition metal atoms surrounded by organic 
ligands. 1 ^* In a SMM crystal, these ligands separate the mag- 
netic core of each molecule from those of its neighbors mak- 
ing the inter-molecule magnetic exchange interactions negli- 
gible. Thus, the magnetic properties of a SMM crystal are 
due to the magnetic moments of the individual SMMs. The 
SMM's large spin and magnetic anisotropy barrier (MAB) 
hamper magnetization reversal below the SMM's blocking 
temperature. Mni2, the first discovered SMM 5 , contains 
four Mn 4+ ions (henceforth referred to as the "inner Mn") 
surrounded by a nonplanar ring of eight Mn 3+ ions (the 
"outer Mn") that interact antiferomagnetically with the inner 
Mn. The large ground state spin S=10 6 and large magnetic 
anisotropy barrier (up to ~ 6.1 meV along the easy axis) 7 
lead to high blocking temperatures (~ 3.5 K) 8 and long relax- 
ation times. 9 For these reasons the Mni2 family has been the 
most studied among SMMs. 

Electron and spin transport through individual molecules 
bridging a pair of electrodes, and through single molecule 
transistors that include a third "gate" electrode, have been 
studied intensively for more than a decade. 10 Because of their 
large magnetic anisotropy barriers and associated stable mag- 
netic moments, single molecule magnets bring a new dimen- 
sion to this field and also raise the possibility of molecular 
magnetic information storage. Therefore the transport prop- 
erties of transistors based on individual SMMs are attracting 
considerable interest at present, both experimentally 11-13 and 
theoretically. 1 1-25 



Transport experiments on individual Mni2-based SMMs 
have demonstrated Coulomb blockade, negative differen- 
tial conductance and magnetic field dependence of the 



transport 



11,12 



On the theoretical side, two different ap- 



proaches have been adopted in studying transport through in- 
dividual SMMs: 

1. Models based on effective spin-Hamiltonians 1 1-20 have 
yielded many important insights into the behavior of these 
systems. They have the advantages of their conceptual sim- 
plicity and of the transparency of the results that they yield. 
However, although in this approach the spin properties of the 
SMM are taken into account, the other aspects of the elec- 
tronic structure of the molecule that play an important role 
in transport (in particular the spatial details of its molecular 
orbitals) are not considered. 

2. Density functional theory (DFT)-based calculations 21 " 26 
have overcome this limitation by treating both the spin and 
spatial electronic structure of a SMM in a unified way. How- 
ever, DFT transport calculations are time consuming. Thus 
in the case of large molecules it has been necessary to re- 
duce the number of atoms (by replacing the SMM with a 
smaller family member 26 or shortening all or most of the 
ligands 21 24 ) in order to make the DFT approach to transport 
calculations practical. This simplification, although reason- 
able, has limitations: Even though the total magnetization of a 
Mn 12 SMM covered with CH 3 ligands has been predicted not 
to change if these ligands are removed 27 or replaced by hydro- 
gen atoms 28 ' 29 , the local features of the electronic structure at 
each Mn site have been predicted to change significantly. 30 
Furthermore, it has been shown experimentally that ligands 
can affect the magnetic properties of the SMM. 8 Also, in or- 
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der to obtain reasonable values of the energy gap between the 
highest occupied molecular orbital (HOMO) and lowest unoc- 
cupied molecular orbital (LUMO) for the SMMs, it has often 
been necessary to introduce into the DFT-based calculations 
an adjustable Hubbard U parameter whose value is not known 
accurately. 22 ' 23 ' 3 1-35 Thus, while DFT + U calculations have 
produced more realistic Mni2 HOMO-LUMO gaps than pure 
DFT calculations, the DFT + U calculations carried out by 
different groups have yielded opposite spin polarizations for 
the Mn 12 LUMO orbital. 22 ' 23 

In this paper we describe and analyze a model, that incorpo- 
rates both the spin and spatial aspects of the electronic struc- 
ture of SMMs. A preliminary account of the model has been 
presented elsewhere. 36 This model describes the fundamen- 
tal properties of Mn 12 SMMs quite well, yielding calculated 
values of the SMM total spin, the magnetic moments of the in- 
ner and outer Mn atoms, the magnetic anisotropy barrier and 
the HOMO-LUMO gap that are consistent with experiment. 
However, it is much simpler than DFT and transport computa- 
tions that are based on it are much less time consuming. Thus 
we have been able to investigate two members of the Mni2 
SMM family including all of their ligands. 

In our previous paper 36 we considered only the neutral 
Mni2-benzoate single molecule magnet. Here in addition to 
discussing the electronic and spintonic structure and transport 
properties of Mni2-benzoate in more detail, we shall consider 
also Mni2-acetate, both its neutral and negatively charged 
species. Mni 2 -acetate and Mni2-benzoate have significantly 
different properties. In particular for Mn 12-acetate the LUMO 
is located on the Mn core of the molecule while for Mni2- 
benzoate we predict the LUMO to be located on the ligands. 
Therefore comparing and contrasting our results for these two 
molecules as we shall do here is of interest. We note that dif- 
ferent theories 22,23 of Mni2-acetate have produced differing 
results regarding its negatively charged species and HOMO- 
LUMO gap. Our theory provides a simple and complete de- 
scription of its properties that is consistent with experiment. 
For neutral and negatively charged Mni2-acetate our tight- 
binding model yields the total spin of the SMM, the spins of 
the individual Mn ions, the magnetic easy axis orientation, 
the size of the magnetic anisotropy barrier and the size of the 
HOMO-LUMO gap consistent with experiment 37 . Therefore, 
the results that we present here for Mn 12-acetate provide a 
useful benchmark for our tight-binding model. 

In the systems considered in the previous theoretical trans- 
port studies of Mni2 SMMs the ligands have simply played 
the role of tunnel barriers and thus the SMMs have always 
been in the Coulomb blockade regime. Our calculations reveal 
that this need not always be the case, with important conse- 
quences for the SMM's transport properties. In particular, for 
a Mni2 SMM covered with 4-(methylthio)benzoate ligands 
we predict that the LUMO and near LUMO orbitals are lo- 
calized not on the Mni2 core of the molecule (as in the SMMs 
considered in previous studies) but on some of the ligands. As 
a consequence we predict that, for some orientations of the 
SMM relative to gold electrodes and positive gate voltages, 
the mechanism of conduction through the molecule is off- 
resonant or resonant tunneling that is not subject to Coulomb 



blockade. The HOMO, however is predicted to be localized 
on the Mn 12 core of the molecule. As a consequence, for neg- 
ative gate voltages transport through the SMM is predicted to 
be in the Coulomb blockade regime. 

In the cases where we predict transport not to be subject to 
Coulomb blockade the magnetic easy axis of the molecule is 
oriented approximately parallel to the direction of the current 
flow through the molecule from one electrode to the other. 
Thus our theory opens the way to resolving for the first time 
an important but previously intractable experimental problem, 
that of determining the direction of the total spin vector of the 
molecule relative to the electrodes (and to the direction of the 
current flow): In the Mni2 SMM transistors that have been 
realized experimentally to date, there has been no control over 
the orientation of the SMM relative to electrodes and how the 
orientation of the easy axis might be identified experimentally 
in these systems has remained an open question. We predict 
that if transport through the molecule is observed to be in the 
Coulomb blockade regime for negative gate bias but not in the 
Coulomb blockade regime for positive gate bias, then the easy 
axis of the molecule is approximately parallel to the direction 
of the current flow. 

Previous theoretical work 24 has suggested that the magni- 
tude of the conductance of a Mni2 SMM molecule may de- 
pend on the orientation of the easy axis relative to the direc- 
tion of current flow. We also find the spin-unresolved current 
through the molecule to depend on the molecule's orientation. 
However, the conductance also depends on the atomic-scale 
details of the bonding between the molecule and electrodes. 
Therefore observation of a qualitative signature (the molecule 
switching into and out of the Coulomb blockade regime de- 
pending on the sign of the gate voltage) as predicted here 
would be a much more convincing way to determine the ori- 
entation of the easy axis. 

It has also been suggested 17 that it may be possible to de- 
termine the orientation of the SMM easy axis experimentally 
by mapping the behavior of transport resonances as functions 
of the bias voltage and direction of a strong external magnetic 
field applied to the system. Detecting instead the presence or 
absence of Coulomb blockade as a function of gate voltage as 
is proposed here will be more straight forward experimentally. 

Our model also predicts strong spin filtering by the SMM 
in the case of a Mni2 molecule in its ground state and strongly 
coupled to gold electrodes. Whether majority spin (up) or mi- 
nority spin (down) electrons are transmitted preferentially de- 
pends on both the sign of the gate voltage and on the particular 
ligands via which the molecule bonds to the electrodes. 

In Section II we introduce our SMM Hamiltonian. Our 
starting point is the extended Hiickel model, a tight binding 
approach that has been used successfully to describe the elec- 
tronic structures of many non-magnetic molecules. 10 We show 
in this Section how the extended Hiickel model can be gener- 
alized so as to describe the magnetism of SMMs, including 
the spin polarization and the effects of spin-orbit coupling. In 
order for the model to describe the SMM magnetic anisotropy 
barrier correctly, we find it to be necessary to include in it 
/nfer-atomic contributions to the spin-orbit coupling operator 
as well as the intra-atomic contributions. We derive analytic 
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tight-binding approximations for both contributions in Section 
II C and the Appendix A. In Section III we present the results 
of our calculations for the magnetic anisotropy barriers and 
spin-resolved densities of states for some Mni2 SMMs, and 
examine how these and the nature of the molecular HOMO 
and LUMO depend on the ligands, the details of the molec- 
ular geometry and the molecular charge. In Section IV A we 
outline the methodology used in our calculations of transport 
in individual SMMs strongly coupled to the electrodes. We 
discuss the relationship between the orientation of the SMM 
easy axis, the SMM molecular orbitals and transport in Sec- 
tion IV B. The lifting of the Coulomb blockade and how this 
may be used in determining the orientation of the magnetic 
easy axis are discussed in Section IV C. We present our pre- 
dictions regarding spin filtering by Mni2 SMMs in Section 
IV D. Our conclusions are presented in Section V. 

II. THE MOLECULAR MODEL HAMILTONIAN 

Our SMM Hamiltonian 16 is a generalization of the semi- 
empirical extended Huckel model Hamiltonian. 10 ' 38,39 It con- 
tains the following terms: 

Here H EH is the extended Huckel Hamiltonian, 10 - 38 - 39 iJ s P in 
gives rise to the spin polarization of the molecule and H so is 
the spin-orbit coupling term. 

A. Extended Huckel model 



B. Spin polarization 

While extended Huckel theory has been used to help 
construct models describing spin transport in non-magnetic 
molecules contacted by magnetic electrodes, 5 1-53 the ex- 
tended Huckel Hamiltonian _ff EH does not by itself take ac- 
count of spin polarization. It also has only one set of parame- 
ters for manganese, whereas in the Mni2 SMMs the inner and 
outer Mn ions have different (+4, +3) oxidation states and 
antiparallel spins. In our model, H spm addresses these issues. 
Its matrix elements (i, s, a\H spm \i', s', a'} = F^°, s , a , be- 
tween valence orbitals i and i' of atoms a and a' with spin s 
and s' are defined as follows: 

H- -, i , = Di a -ii a i (s\n ■ Sis > 

{dinner if a is an inner Mn and i is a d- valence orbital 
Pouter if ot is an outer Mn and i is a rf-valence orbital 
otherwise 

(2) 

Here ^i nner and ^4. outer are fitting parameters, n is a unit vector 
aligned with the magnetic moment of the SMM, and S is the 
one-electron spin operator. In the Mni2 ground state the total 
spin of each inner Mn (Sinner = — §■) is antiparallel to the 
total spin of each outer Mn (S ou tei = 2). 2 Also the total spin 
of the Mn 12 SMM is parallel to the spins of the outer Mn ions 
and antiparallel to the spins of the inner Mn ions. 2 Therefore 
Annei- and Aouter have opposite signs (A inner > 0, A outa - < 0). 
Numerical values of these parameters are given in Sec. Ill A 
and Sec. IIIB. 



In extended Huckel theory, a small set of Slater type atomic 
valence orbitals \^i a ) is chosen as the basis. Here \^i a ) 



is the i atomic orbital of the a atom. The extended 
Huckel Hamiltonian's diagonal elements i a \H Eil \^ i a ) = 
Hfa-ia = £ ia are the experimentally determined negative va- 
lence orbital ionization energies £j a . The non-diagonal ma- 
trix elements are given by H™ ila , = D ia . i ' a 'K £io:+ ^ i ' a ' , 
where Di a: i> a r = (^i a \^i'a') are the orbital overlaps and K 
is an empirical factor, chosen for consistency with experimen- 
tal molecular electronic structure data. In our calculations, 39 
K = 1.75 + A 2 ia;i , a , - 0.75Af Q .,, Q , where A ia ^ a , = 
(e ia - £i' <*')/(£ ) as was proposed in Ref. 38. 

In recent years, extended Huckel theory has been used 
to understand the transport properties of a variety of non- 
magnetic molecular systems: Calculations based on extended 
Huckel theory have yielded tunneling conductances 40-43 and 
inelastic tunneling intensities 44 in agreement with experiment 
for molecules thiol bonded to gold electrodes. They have 
also explained transport phenomena observed in STM exper- 
iments on molecular arrays on silicon 45 ^ 7 as well as elec- 
troluminescence data 48 , current-voltage characteristics 48 and 
STM images 49 of molecules on complex substrates. Re- 
cently, extended Huckel theory has also been used success- 
fully to model the electronic structures of adsorbates cova- 
lently bonded to graphene and electron transport in graphene 
nanoribbons with such adsorbates. 50 



C. spin-orbit coupling 

Spin-orbit coupling is also not included in extended Huckel 
theory. However, it plays a very important role in SMM 
physics since it is responsible for the magnetic anisotropy of 
these systems. We find that in order to describe the SMM 
magnetic anisotropy correctly within a tight binding model 
it is necessary to consider the /nfer-atomic contributions to 
the spin-orbit Hamiltonian operator H so (perturbation of the 
atomic orbitals by the core potentials of other atoms) as well 
as ;rafra-atomic contributions (perturbation of the atomic or- 
bitals by the core potential of the same atom). We derive 
appropriate analytic expressions for these matrix elements in 
this paper. Our results generalize extended Huckel theory to 
include the effects of spin-orbit coupling. Spin-orbit coupling 
in SMMs has previously been treated within the framework 
of density functional theory. 28 - 54 ~ 56 The present tight-binding 
theory has the advantages that it is analytic and relatively 
straight forward. 

The spin-orbit Hamiltonian is given by 57 



H 



so 



-a ■ W(r) x p 



(3) 



(2mc) 2 

where p is the momentum operator, V(r) is the electron 
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Coulomb potential energy, cr = (<r x ,a y , a z ) and <r x , a y and 
<j z are the Pauli spin matrices. 

We approximate V(r) by a sum of atomic electron poten- 
tial energies V(r) ~ S a K«(r — r a ) where r a is the po- 
sition of a th atomic nucleus. Since the spin-orbit coupling 
arises mainly from the atomic cores where the potential en- 
ergy V a (r — r a ) is approximately spherically symmetric we 



approximate VV^(r — r a ) 
becomes 



dV a (\r-r a \) 
|r— r a | d\r~ r n | ' 



Then Eq. 3 



H 



so 



E 



i 



dV a {\r-r a \) 



2m 2 c 2 |r — r Q | d[r — r a | 



S • L Q (4) 



where S = \ ficr is the electron spin angular momentum oper- 
ator, L Q = (r — r Q ) x p is the atomic orbital angular momen- 
tum operator with respect to the position of the a th nucleus 
and the sum is over all atoms a. The matrix elements of H so 
between valence orbitals i and %' of atoms (3 and /?' with spin 
s and s' are evaluated in Appendix A and we find 



^IH^'s'P') ~ E&j.fi&w + (1 - Sp, 

3 

+ [Di'P'tipEjXs-p]*) 



(5) 



Here the first term on the right-hand side is the intra- 
atomic contribution and the remaining terms are the inter- 
atomic contribution. An explicit expression for E™^f s ,.p, the 
spin-orbit coupling Hamiltonian matrix element due to intra- 
atomic spin-orbit coupling, is Eq. A2. 



III. PROPERTIES OF ISOLATED Mm 2 MOLECULES 

In this Section we present the predictions of the model in- 
troduced in Section II for isolated SMM molecules that are 
based on solutions of the tight binding Schrodinger equation 



H 



SMM 



(6) 



for the molecular orbitals, \<j>i), and their energies, £",. 

It should be noted that in the present work the values of 
the parameters H™ ila , and -Di a ,j' a ' that enter the extended 
Hiickel model, the spin polarization Hamiltonian, Eq. (2), and 
the spin-orbit coupling Hamiltonian, Eq. (A5), were adopted 
without modification from Refs. 38 and 39. 

Experimental estimates 60 of the spin-orbit coupling con- 
stant for manganese, eMn,h (see the discussion following Eq. 
A2 in Appendix A), have been in the range 0.023-0.051 eV, 
while theoretical estimates 59 have been in the range 0.038- 
0.055 eV. In this paper for manganese atoms we use the value 
tyin.h =0.036 eV, which is consistent with the experimental 
and theoretical values. Our calculations show that the spin- 
orbit coupling constants of the other atoms in molecule do not 
effect the SMMs' properties significantly. However in the nu- 
merical results presented below the spin-orbit coupling due to 
the O and Au atoms in the extended molecule as well as that 



due to the Mn atoms is included. Thus the only free parame- 
ters in the present model are A mn ei and ^4. outer of Eq. (2) that 
control the spin polarizations of the Mn atoms. 



A. Magnetic anisotropy barrier 

The simplest SMM that we shall discuss here is 
Mni 2 Oi2(02CR)i 6 (H 2 0)4 with R = CH 3 , henceforth referred 
to as Mni2-Ac. The molecular geometry of Mni2-Ac used in 
the present calculations was obtained from the experimentally 
measured 61 geometry of Mni 2 Oi2(02CR)i 6 (H 2 0)4 for R = 
CHCHCH 3 by substituting CH 3 for R. 

For the Mni2-Ac molecule we chose the values of the mag- 
netic parameters in Eq. 2 to be Anna- =3.0 eV and Amer = 
-2.6 eV. For these values, we find the calculated ground state 
total spin of the neutral Mni2-Ac to be S to tai = 10 in agree- 
ment with the experimental value. 62 We find the calculated 
local magnetic moments associated with inner and outer man- 
ganese ions to be -3.18/Ub and 3. 86/13, respectively. Accord- 
ing to Hund's rules, these magnetic moments (that correspond 
to spin — I and +2, respectively) are in a good agreement with 
the Mn 4+ (Mn 3+ ) oxidation states of the inner (outer) man- 
ganese ions. 5 They are also comparable with the results of 
other theoretical calculations. 27 ' 28 

We estimate the magnetic anisotropy energy of the SMM 
from the total ground state energy expression 



(7) 



where Ei and |</>j) are molecular orbital (MO) eigenenergies 
and eigenstates as in Eq. (6) and the summations are over 
all occupied MOs. Since H spm represents electron-electron 
interactions at a mean field level, the second summation on the 
right-hand side of Eq. (7) is required to avoid double counting 
the corresponding interaction energy. 

In order to evaluate the magnetic anisotropy barrier (MAB) 
of the SMM, the total energy of the molecule, Eq.(7), has been 
calculated as a function of the projection of the total spin on 
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Figure 1. (Color online) Calculated energy of the neutral Mni2-Ac 
molecule with ground state spin Stotai = 10 versus the total spin 
projected on the easy axis. The arrow shows the magnetic anisotropy 
barrier. 
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the easy axis. As shown in Fig. 1, the calculated MAB of neu- 
tral Mni2-Ac is 6.25 meV which is close to its experimental 
value of 6.1 meV 7 . Our calculations show the contributions of 
/nfer-atomic and intra-atomic spin-orbit coupling to the MAB 
to be of the same order of magnitude. The results presented 
above demonstrate that our tight binding model is able to pro- 
vide a very good description of the experimentally observed 
magnetic properties of the Mni2-Ac SMM. 

B. Role of the ligands 

As has already been mentioned, in DFT-based calculations 
for SMMs, most or all of the ligands have usually been re- 
placed by hydrogen atoms and thus ligand effects have not 
been examined fully. 23 26 In order to investigate the ligands' 
influence on SMM properties, we considered another member 
of the Mn 12 family, Mn 12 Oi2(02CC 6 H5) 16 (H 2 0)4 (referred 
to as Mni2-Ph), that has larger (benzoate) ligands than Mni2- 
Ac. 36 In the present work the geometry of Mni2-Ph has been 
taken from experimental data. 63 In our modeling of Mni2-Ph 
we chose the values of the magnetic parameters in Eq. 2 to 
be Annei- =3.0 eV and A>uter =-3.5 eV. For these parameter 
values we find the calculated ground state's total spin to be 
S total = 10 and the calculated MAB to be 2.50 meV (see Fig. 
4 (a)), values that are consistent with experiment. 63 ' 64 

To compare the effects of the different ligands for Mni2, 



we have plotted the calculated densities of states projected on 
the oxygen, outer Mn, inner Mn atoms, and carbon atoms in 
Figs. 2 and 3 for Mni2-Ac and Mni2-Ph, respectively. Up 
to the HOMO energy Fig. 2 and Fig. 3 are similar for both 
the inner and outer Mn. In both cases outer Mn are filled 
mostly with spin up electrons (parallel to the total spin), con- 
sistent with having spin S ou ter Mn = 2, the inner Mn are filled 
with spin down electrons (anti-parallel to the total spin) in 
agreement with Si nne r-Mn = — § > an d the carbon atoms are not 
strongly spin polarized. As is seen in in Figs. 2 and 3, the 
HOMO level in both cases is predicted to be spin up (major- 
ity spin) and located mainly on the outer Mn and on the O 
atoms, consistent with the results of recent generalized gra- 
dient approximation (GGA) and local density approximation 
(LDA)+U density functional calculations for Mni2-Ac that in- 
clude the complete set of CH3 ligands. 23 

The calculated energy gaps between the HOMO and 
LUMO (lowest unoccupied molecular orbital) for both the 
Mn 12 -Ac and Mn 12 -Ph molecules are on the order of leV 
which is common for Mn 12 family. 65 It is worth mentioning 
that in DFT calculations, it is typical that LDA and GGA un- 
derestimate the HOMO-LUMO gaps in transition metal oxide 
systems 66 and that this has been corrected by using the on- 
site LDA+U correction where the Hubbard U is treated as 
an adjustable parameter. For instance the calculated HOMO- 
LUMO gap for Mni 2 -Ac increased from 0.35 eV in the GGA 
to 1.6 eV in LDA+U in some theoretical studies. 23 
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Figure 2. (Color online) Projected density of states for majority spin 
(up) and minority spin (down) electrons on the oxygen, inner Mn, 
outer Mn and carbon atoms for Mni2-Ac. 



Figure 3. (Color online) Projected density of states for majority spin 
(up) and minority spin (down) electrons on the oxygen, inner Mn, 
outer Mn and carbon atoms for Mni2-Ph. 
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The LUMO can give an indication as to where the 
added electron may be located in a negatively charged 
molecule and of the negatively charged molecule's prop- 
erties. Several experiments have been reported in which 
reduced anionic products of Mni2 have been synthesized 
in the presence of the PPlu (PPlu = P(C6Hs)4) molecule 
which acts as an electron donor and their properties 
have been studied. 37,68 ' 69 Among the studied molecules are 
Mni 2 Oi 2 (0 2 CEt)i 6 (H 2 0) 4 37 , where Et = CH 2 CH 3 , (hence- 
forth Mn 12 -Et) and Mn 12 -Ph 68 . To our knowledge there has 
not been any experimental study of a reduced product of 
Mni 2 -Ac. Since Mni 2 -Et has only one more methylene in 
its ligands than Mni 2 -Ac, it is reasonable to assume that they 
have similar properties. 

As is seen in Fig. 2 and 3, our model predicts the LUMOs of 
both Mni 2 -Ac and Mni 2 -Ph to be spin down electron states, 
although for Mni 2 -Ph there are also spin up states that are 
close in energy to the LUMO (see also Fig. 10 ). Based on 
this the calculated total spins of the negatively charged Mni 2 - 
Ac and Mni 2 -Ph molecules are both predicted by the present 
model to be S total = 9^, in agreement with experiments. 37 ' 68 
Although in both cases the LUMOs are spin down states, the 
predicted locations of the LUMOs for the two molecules are 
different: In Mni 2 -Ac, an added electron is located primarily 
on the outer Mn's as is seen in Fig. 2, whereas in Mni 2 -Ph, 
the location of an added electron is primarily on carbon atoms, 
see Fig. 3. 

It has been found experimentally 37 that an added electron 
on Mni 2 -Et is localized on an outer manganese and the extra 
electron produces a trapped-valence Mn u Mn!- u Mn4 V system. 
Our calculation yielded a similar result for Mni 2 -Ac, where 
the LUMO is localized primarily on outer manganese atoms. 
For Mni 2 -Ac, the present model predicts the energy gaps of 
majority and minority spin to be 0.83 and 3.1 eV, respectively. 
These values are in good agreement with corresponding GGA 
estimates 23 ' 28 of 0.4 and 2.1 eV for majority and minority 
spin energy gapes, considering that GGA underestimates the 
HOMO-LUMO gaps in transition metal oxide systems. 66 For 
Mni 2 -Ac the GGA also yields a spin up LUMO primarily on 
the O and outer Mn while LDA+U yields a spin down LUMO 
primarily on the inner Mn and O. 23 Interestingly, the result 
obtained from the GGA 23 for the location of the Mn 12 LUMO 
appears to be in better agreement with both the prediction of 
our model and the above-mentioned experiment 37 than that 
obtained from LDA+U. 2j Therefore, the qualitative features 
of the HOMO and LUMO of Mn i2 -Ac in the present model 
agree with experiment on the related Mni 2 -Et molecule, and 
more consistently than of the DFT-based calculations in the 
current literature. 

For negatively charged Mni 2 -Ph it was found to be too 
difficult to determine the valence of each manganese ion 
experimentally. 68 Therefore experimental information regard- 
ing where the extra electron is localized in this molecule is not 
available in this time. Our calculation predicts that in Mni 2 - 
Ph the added electron is localized on the benzoate ligands (in- 
cluding carbon and oxygen atoms) instead of the manganese 
atoms. To our knowledge DFT-based calculations of the elec- 
tronic structure of the Mni 2 -Ph molecule with the complete 



set of ligands are not available in the literature at this time. 



C. Role of geometry 

To investigate the properties of neutral Mni 2 SMMs, we 
have used the available experimental geometries of the neu- 
tral Mni 2 molecules. For negatively charged molecules, 
as a first approximation, we have assumed in the preced- 
ing Section (IIIB) that the added electron locates where 
LUMO of the neutral molecule is located, without changing 
the molecule's geometry. We shall call this the "unrelaxed 
charged molecule". By making this idealization we were able 
to calculate the spin and magnetic anisotropy barrier (MAB) 
of the negatively charged SMMs. 

As has already been mentioned, there have been several 
experiments studying reduced products of Mni 2 . 37,69 There- 
fore as a better approximation, we have also carried out cal- 
culations for the Mni 2 -Ph molecule using the experimen- 
tally determined geometry of the negatively charged Mn 12 -Ph 
molecule 68 . We will call this geometry the "realistic charged 
molecule geometry." 

Fig. 4 (a,b) show the calculated MABs of the neutral and 
unrelaxed negatively charged Mni 2 -Ph. Even though, accord- 
ing to the experimental data, 68 the geometry of the negatively 
charged Mni 2 -Ph differs from that of the neutral molecule, as 
is seen in Fig. 4 (b) and (c) our calculations did not show a 
large change in the MAB due to this difference. This insensi- 
tivity of the MAB to the charge state is in agreement with ex- 
periment: The experimental values of the lower-temperature 
MAB for the neutral 64 and negatively charged molecule 70 are 
3.3 meV and 2.41 meV, respectively. So, experimentally, the 
MAB changes by 27% if an extra electron is added. 

Our model offers a physical reason for the insensitivity of 
the MAB to the charge state, as shall be explained next: As 
is seen in Fig. 3 the LUMO of Mni 2 -Ph is predicted to be lo- 
cated mainly on carbon atoms and not on the manganese. Our 
calculations yielded a similar result for the realistic charged 
molecule geometry. On the other hand the main source of the 
MAB of Mn i2 is the Jahn-Teller distortion of the Mn+ 3 , see 
Ref. 12. Therefore even though the geometry has changed 
due to addition of an extra electron, because the extra elec- 
tron does not change the oxidation state of the Mn atoms the 
Jahn-Teller distortion does not change significantly and conse- 
quently the MAB does not change by much when an electron 
is added to the molecule. 



IV. TRANSPORT 

A. Theory of coherent transport in the strong coupling regime 

In electronic transport experiments the SMM is connected 
to source and drain leads. Thus, the Hamiltonian for the entire 
system may be written as: 

H = H EM + H R + H h + W R + W h (8) 
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Figure 4. (Color online) Calculated energies of Mni2-Ph versus total 
spin projected onto the easy axis for (a) neutral molecule with spin 
ground state S tota l — 10 , (b) unrelaxed (explained in text) negatively 
charged molecule with S tota i = 9|, and (c) realistic (explained in 
text) negatively charged molecule with S tot ai = 9§- Arrows show 
the magnetic anisotropy barriers. 



Here the system has been divided into a semi-infinite right 
lead H R , a semi-infinite left lead H h and an extended 
molecule H EM . The extended molecule contains the central 
SMM as well as a cluster of nearby atoms (assumed here to 
be gold) belonging to each of the two contacts. We have car- 
ried out systematic calculations to determine an adequate gold 
cluster size for which the transport properties are not affected 
by the size of the clusters. The results that we present below 
were obtained for extended molecules with 3 1 gold atoms in 
each cluster. There cluster sizes were found to be adequate in 
the above sense. V^ R(L) is the coupling Hamiltonian between 
the infinite right(left) lead and extended molecule. 

To evaluate H EM we have used Eq. (1) applied to the ex- 
tended molecule. The right and left leads were modeled as a 
large number of semi infinite one-dimensional ideal channels 
that represent macroscopic electron reservoirs, as in previous 
studies of electron and spin transport through single molecules 
with gold contacts. 43 ' 44 ' 7 1-74 Nine semi-infinite ideal atomic 
chains attached to each of the eighteen gold atoms of each 
cluster that are furthest from the SMM represent the source 
and drain leads in the transport results that we shall present 
here. 

2 

According to the Landauer formula, g = ^-T(Ep), the con- 
ductance is proportional to the total transmission probability 
T at the Fermi energy Ep. For non-zero temperatures and ap- 
plied bias voltages V, Landauer theory yields 10 



I(V) 



dET(E,V)(f(E-f, a )-f(E-f, d )) (9) 



where / is current, f(e) is the Fermi function, /i s (/^) is the 
electrochemical potential of the source(drain) electrode. The 
total transmission probability can be written as 



T(E,V) = J2 U f- 

IJSS' 



js';is 



(10) 



E is the incident electron's energy, ij S ' : ; s is the transmission 
amplitude from i th electronic channel of left lead with spin 
s and velocity Vi S to j th electronic channel of the right lead 
with spin s' and velocity Vj S /. 



We find the transmission amplitudes t by solving the 
Lippmann-Schwinger equation 



I*") = \®o) + G (E)W\$ a ) 



(11) 



where |$q ) is the eigenstate of the a decoupled semi infinite 
one dimensional lead, Go(E) is the Green's function of the 
decoupled system, W is coupling matrix between extended 
molecule and leads, and \*f? a ) is the scattering eigenstate of 
the coupled system associated with the incident electron state 
$q). Since the basis set used in extended Huckel theory is 
non orthogonal, we apply the orthogonalization procedure de- 
scribed in Ref. 75 in these calculations. 

The Landauer theory described above neglects charging 
effects that can give rise to Coulomb blockade in a single 
molecule magnet that is very weakly coupled to all of the 
leads that carry electrons to and from the molecule. In some 
of the cases to be considered below this coupling is found to 
be very strong so that Coulomb blockade should not occur 
and the above Landauer theory is expected to be valid. How- 
ever, even in the cases where the coupling is weak, calcula- 
tions based on the Landauer theory provide useful qualitative 
insights into the molecular states that participate in transport, 
the degree of broadening of these levels due to their coupling 
with the leads and hence whether or not Coulomb blockade 
may be expected to occur. 



B. Interplay between the orientation of the magnetic easy axis 
relative to the leads and electronic structure and transport. 

In the experiments in which transport through a SMM has 
been measured 1 1-13 there has been no control over the orien- 
tation of the SMM relative to the electrodes. In particular, the 
orientation of the SMM easy axis has not been controllable. 
While experimental control over the orientation of the SMM 
easy axis has not been achieved, it is of interest to consider 
whether it is possible to experimentally determine the orien- 
tation of the easy axis relative to the electrodes that is realized 
in individual samples by making appropriate measurements. 

We have investigated this by calculating the transport prop- 
erties of SMMs molecules bridging gold electrodes in a va- 
riety of bonding configurations that correspond to different 
orientations of the easy axis relative to the electrodes. For 
our studies of coherent quantum transport through SMMs, 
we have chosen 4-(methylthio) benzoate Mn 12 henceforth re- 
ferred to as Mni2-Ph-Th. That is, all of the benzoate ligands 
are terminated with methylthio (SCH3) groups (see Fig. 5) 
except for two ligands that are terminated with sulfur atoms 
(with their methyl removed) that bond chemically to the gold 
leads. 

Fig. 6 shows the projected density of states of Mni2-Ph- 
Th on the sulfur, oxygen, outer Mn, inner Mn, and carbon 
atoms. The Mni2-Ph-Th projected density of states is similar 
to that of Mni2-Ph in Fig. 3. However there are noticeable 
differences such as the appearance of the sulfur and carbon 
density of states features below the HOMO in Fig. 6, as well 
as the wave functions of the LUMO and MOs near the LUMO 
in energy penetrating from the carbon to the sulfur atoms. 
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If the Mni2-Ph-Th molecule shown in Fig. 5 is oriented 
so that its magnetic easy axis is along the z-axis then the Mn 
atoms are located near the x-y plane. With the coordinate axes 
chosen in this way, we refer to the ligands that are close to the 
x-y plane as "in-plane ligands" and to the others as "out-of- 
plane ligands". We considered the bonding configurations la- 
belled IPBi - IPB6 in Fig. 5 (a) in which the molecule bonds 
to the leads via different pairs of in-plane ligands and the con- 
figurations OPBi - OPB4 in Fig. 5 (b) where the molecule 
bonds to the leads via pairs of out-of-plane ligands. 

Fig. 7 (Fig. 8) shows the total electron transmission prob- 
ability through the molecule as a function of energy, E, rela- 
tive to Fermi energy of electrodes 67 , Ep, for the different IPB 
(OPB) configurations. The details of the calculated transmis- 
sion spectra for the in-plane-bonding configurations shown in 
Fig. 7 depend on the ligands involved in the bonding. How- 
ever, in all cases the transmission probabilities decay very 
rapidly when the energy falls below that of the LUMO or 
rises above that of the HOMO. On the other hand, this is not 
the case for some of the out-of -plane-bonding configurations 
as can be seen in Fig. 8: In particular, for OPBi and OPB2 
the transmission peaks near the LUMO show pronounced tails 
that are at least two orders of magnitude stronger than those 
for the in-plane bonding configurations in Fig. 7 and extend 
well into the HOMO-LUMO gap. This difference is mani- 
fested in the current carried by the molecule for similar (low 
and moderate) bias voltages applied across the molecule: As 
can be seen in Fig. 9, the calculated currents for the OPBi 
and OPB2 configurations at low and moderate bias are two or- 



1PB 



IPB, 




(b) 



OPB 4 



IPB, IPB fi 



IPB, 




IPB[ OPB 
IPB, 



IPB,, 



OPB, OPB, 



Figure 5. (Color online) Two views of the Mni2-Ph-Th SMM whose 
ligands are terminated with methylthio (SCH3) groups. The Mn 
atoms are located near the x-y plane and the magnetic easy axis is 
aligned with the z-axis. We refer to the ligands that are close to the 
x-y plane as "in-plane ligands" and to the others as "out-of-plane 
ligands". The molecule is assumed to be located between two gold 
leads (not shown) and to bond to each of the gold leads via the sul- 
fur atom of a single ligand, the methyl group having been removed 
from that sulfur atom. The pairs of open circles with the same color 
indicate pairs of sulfur atoms that are attached to the gold leads in 
the different bonding configurations that we consider, (a) The red, 
blue, green, orange, purple, and brown circles labelled IPBi indi- 
cate the pair of sulfur atoms bonding to the gold in the i bonding 
configuration involving in-plane ligands. (b) The red, green, purple, 
and orange circles labelled OPBi indicate the pair of sulfur atoms 
bonding to the gold in the i th bonding configuration involving out- 
of-plane ligands. Atoms are color labeled: manganese (red), carbon 
(grey), sulfur (yellow), oxygen (blue), and hydrogen (white). 



ders of magnitude larger than for the IPB configurations. On 
the other hand, the currents for the OPB 3 and OPB 4 config- 
urations (for which the tails of the transmission peaks below 
the LUMO in Fig. 8 are much weaker than for OPBi and 
OPB2) are similar in magnitude to those for some of the IPB 
configurations. These differences between the transport prop- 
erties of the SMM in the various bonding geometries may be 
understood physically as follows: 

The LUMO and the other molecular orbitals that lie closest 
to it in energy are located mainly on the out-of-plane ligands 
and have much weaker weight on the in-plane ligands (see 
orbitals (a) and (b) in Fig. 10). This favors a strong hybridiza- 
tion of some of the near-LUMO molecular orbitals to gold 
electrodes that bond to the molecule via out of plane ligands. 
In cases where this strong hybridization occurs, off-resonant 
transport at energies below the LUMO through molecules 
bound to the electrodes via out-of-plane ligands is expected 
to be much stronger than that through molecules bound via 
in-plane ligands. This is consistent with the much larger cur- 
rents seen in Fig. 9 for OPBi and OPB 2 configurations than 
for all of the IPB configurations. Note that the molecular or- 
bitals with larger weight on the in-plane ligands are located 
further in energy from the LUMO. They lie at least 0.18 eV 
above the LUMO as is shown in Fig. 10 (c). Thus off-resonant 
transport below the LUMO via these orbitals even for gold 
electrodes bonding to the in-plane ligands is expected to be 
relatively weak as is seen in Fig. 9. However, other consider- 
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Figure 6. (Color online) Projected density of states for majority spin 
(up) and minority spin (down) on the sulfur, oxygen, inner Mn, outer 
Mn, and carbon atoms for Mni2-Ph-Th. 
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Figure 7. (Color online) Total transmission probability at zero bias 
and gate voltage as a function of energy, E, relative to Fermi energy, 
Ef, of gold for the different in-plane bonding geometries between 
the molecule and leads that are identified in Fig. 5 (a). The inset is 
a schematic of the SMM and electrodes. The magnetic core of the 
molecule that contains the Mn atoms is pink, the arrow indicates the 
magnetic easy axis and the electrodes are yellow. 



ations also play a very important role in determining how well 
the molecule conducts, as will be discussed next. 

As can be seen in in Fig. 5, every ligand is attached via 
oxygen atoms to two Mn atoms. We will refer to these as the 
"terminating Mn atoms" of the ligand. We find the magnitude 
of the current through the molecule at low and moderate bias 
to be strongly affected by the distances between the terminat- 
ing Mn atoms belonging to the two ligands that connect the 
molecule to the electrodes and also by the numbers of oxygen 
atoms that bond simultaneously to the terminating Mn atoms 
of those two ligands. These oxygen atoms can be regarded as 
electronic pathways connecting the two ligands. 

For the OPBi geometry two of the terminating Mn atoms of 
the different ligands are both inner Mn separated by 2.819 A. 
In addition to their small separation, these two Mn atoms are 
connected to each other via two oxygen atoms, i.e., two elec- 
tronic pathways. For OPB2 the two ligands that bond to the 
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Figure 8. (Color online) Total transmission probability at zero bias 
and gate voltage as a function of energy, E, relative to Fermi energy, 
Ef, of gold for the different out-of-plane bonding geometries be- 
tween the molecule and leads that are identified in Fig. 5 (b). The 
inset is a schematic of the SMM and electrodes. The magnetic core 
of the molecule that contains the Mn atoms is pink, the arrow indi- 
cates the magnetic easy axis and the electrodes are yellow. 
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Figure 9. (Color Online) Calculated current (from Eq. 9) as a func- 
tion of bias voltage at zero gate voltage and temperature for different 
bonding configurations of the SMM relative to electrodes identified 
in Fig. 5. 



electrodes share the same (outer) terminating Mn atom. The 
distance between the other two terminating Mn atoms (one of 
them is an inner Mn) is also small (2.765 A), however they are 
connected to each other through only one oxygen atom (one 
pathway). For OPB3 all terminating Mn atoms are outer Mn. 
The smallest distance between two terminating Mn atoms of 
the different ligands for OPB 3 is 3.393 A and these Mn atoms 
are only connected to each other through one oxygen (one 
pathway) resulting in the much weaker off resonant transport 
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seen for OPB 3 than for OPBi and OPB 2 in Fig. 9. For OPB 4 
one of the terminating Mn atoms is an inner Mn. The smallest 
distance between two terminating Mn atoms of the different 
ligands for OPB4 is 5.44 A and there is no oxygen atom that 
bonds simultaneously to these terminating Mn atoms (no di- 
rect pathway) resulting in the much weaker off resonant trans- 
port seen for OPB 4 than than for OPBi and OPB 2 in Fig. 9. 

For the in-plane-bonding the terminating Mn atoms are al- 
ways outer Mn atoms. Among them, IPB5 and 1PB$ have 
the smallest distances between their terminating Mn atoms, 
those distances being similar to that for OPB 3 . In addition, for 
OPB3 the terminating Mn are also all outer Mn atoms. These 
resemblances result in the similar currents for OPB 3 and IPB 5 
and IPB 6 (see Fig. 9). 

In the transport calculations presented above, we consid- 
ered SMMs connected to each gold electrode through one 
sulfur-gold atomic junction with a sulfur atom at a hollow 
site (sulfur atom located at the same distance from three gold 
atoms on a [1 1 1] gold surface). In the case of top site bonding 
(sulfur located at the top of a gold atom) the different bond- 
ing geometries have qualitatively similar behavior although 
the currents are two orders of magnitude smaller than for the 
hollow site bonding. We have also found currents enhanced 
by a factor of approximately two if two out-of-plane ligands 
instead of one bond via sulfur atoms to each electrode 36 . 



C. Coulomb blockade versus coherent resonant tunneling and 
the orientation of the magnetic easy axis 

As can be seen in Fig. 10 the LUMO of Mni2-Ph-Th is 
located mainly on some of the out of plane ligands of the 
molecule as are the other molecular orbitals that are closest in 
energy to the LUMO. By contrast the HOMO seen in Fig. 1 1 
is located on the central core of the molecule. Thus if the gold 
electrodes bond chemically to sulfur atoms of the out-of-plane 



ligands on which the LUMO or/and other MOs close in en- 
ergy to the LUMO reside then these MOs may hybridize with 
gold contacts and therefore couple strongly electronically to 
the gold electrodes. We indeed find this strong hybridization 
and strong coupling to the electrodes to occur for some of 
the molecular orbitals that lie close in energy to the LUMO if 
the molecule bonds to the electrodes in the out-of-plane bond- 
ing geometries OPBi and OPB2 for which the molecule is 
the most conductive, as is seen in Fig. 9. Therefore we pre- 
dict resonant and off -resonant transport via these MOs not to 
be subject to Coulomb blockade for those bonding geome- 
tries. We note that this is a novel prediction of the present 
model 36 as in all of the Mni2 SMM systems in previous the- 
oretical transport studies neither the HOMO nor the LUMO 
had significant presence on any of the ligands. Thus the lig- 
ands behaved as tunnel barriers, always resulting in Coulomb 
blockade of transport. 

On the other hand, since in Fig. 1 1 the HOMO and orbitals 
close to it in energy have very little presence on the ligands, 
we predict that the ligands of Mni2-Ph-Th should behave as 
tunnel barriers in the usual way for transport via the HOMO 
and nearby orbitals and therefore transport via these orbitals 
is predicted to be in the Coulomb blockade regime regard- 
less of how the ligands bond to the gold contacts. Further- 
more, if the molecule is oriented relative to the electrodes in 
such a way that only in-plane ligands bond to the electrodes 
then transport via the LUMO and nearby orbitals is also pre- 
dicted to be in the Coulomb blockade regime. Note, however, 
that for some possible bonding configurations between out-of- 
plane ligands and the gold electrodes, hybridization between 
the electrodes and LUMO and near LUMO orbitals need not 
be strong and/or tunneling through the core of the molecule 
may be weak. In such cases Coulomb blockade of transport 
via these orbitals may still occur or transport at energies below 
the LUMO, while coherent, may still be very weak. 

An important consequence of the above predictions is 
that if in a Mni2-Ph-Th transistor transport via the HOMO 
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Figure 10. (Color Online) Projected density of states(PDOS) on car- 
bon atoms for energies close to the LUMO and related wave function 
isosurfaces for Mni2-Ph-Th. The energies of the states with wave 
function isosurfaces (a), (b), and (c) are indicated by arrows in the 
PDOS plot. The LUMO+5, LUMO+14, and LUMO+17 molecular 
orbitals are marked in the PDOS plot with filled red circles. 
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Figure 11. (Color Online) Projected density of states(PDOS) on 
outer Mn atoms in an energy range close to the HOMO and related 
wave function isosurfaces for Mni2-Ph-Th. The energies of the states 
with wave function isosurfaces (a), (b) and (c) are indicated by ar- 
rows in the PDOS plot. 
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is observed to be in the Coulomb blockade regime (i.e., 
Coulomb blockade is observed for negative gate voltages) 
but transport via the LUMO is observed not to be subject to 
Coulomb blockade (i.e., Coulomb blockade of the LUMO is 
not observed for positive gate voltages), then the Mni2-Ph- 
Th molecule is oriented in such a way that at least one out of 
plane ligand is bonded to an electrode. 

Since the easy axis of the Mni2 SMM is roughly aligned 
with the out-of -plane ligands, this means that it should be pos- 
sible to identify specific experimental realizations of Mni2- 
Ph-Th transistors in which the molecular easy axis (and hence 
the molecular magnetic moment) is approximately aligned 
with the direction of current flow through the molecule from 
electrode to electrode simply by observing the presence or ab- 
sence of Coulomb blockade at the thresholds of conduction 
for positive and negative gate voltages. 

It should be realized, however, that different molecular or- 
bitals that are close in energy to the LUMO may be located 
primarily on different out-of-plane ligands. Not all of the lig- 
ands will bond to an electrode in any particular experimen- 
tal realization of a Mni2-Ph-Th transistor. Transport via the 
molecular orbitals that are close in energy to the LUMO but 
are located on ligands that have not bonded to an electrode 
will be subject to Coulomb blockade and its onset will there- 
fore occur at higher positive gate voltages. 



D. Spin filtering 

Previous theoretical work has revealed that Mn 12 molecules 
should act as spin filters even when connected to nonferro- 
magnetic gold electrodes. 16 - 21 ~ 23 In DFT-based studies at the 
level of the generalized gradient approximation (GGA), 22 "" 24 
the Fermi energy of the gold electrodes was found to be be- 
tween the molecular HOMO and LUMO levels which were 
both identified as majority spin (up) states. The minority spin 
(down) MOs were found to be well separated in energy from 
the majority spin (up) MOs and also from the Fermi energy. 
On the other hand, DFT-based calculations at the level of LDA 
+ U yielded a majority spin HOMO but either a majority 22 or 
minority 23 spin LUMO. However, in each case only major- 
ity spin electrons were predicted to be transmitted through the 
Mn i2 SMM at low bias voltages. 22-24 

As was discussed in Section III, the present model predicts 
the HOMO to be a majority spin (up) state and the LUMO to 
be a minority spin (down) state for Mni2-Ac, Mni2-Ph, and 
Mni2-Ph-Th, consistent with experiments. 37 ' 68 This result is 
also qualitatively consistent with one of the previous LDA + 
U theoretical studies that predicted a majority spin HOMO 
and minority spin LUMO for Mni2. 23 However, in the present 
model the gold Fermi level is located near the center of the 
molecular HOMO-LUMO gap 67 for Mm 2 -Ph and Mm 2 -Ph- 
Th. Thus as will be seen below spin-filtering with either ma- 
jority or minority spin electrons being transmitted is predicted 
to be possible for these molecules, depending on the applied 
gate voltage and orientation of the molecule. 

In Fig. 12 we present the results of our spin-resolved trans- 
port calculations for a particular bonding geometry in which 



the two in-plane ligands labelled IPBi in Fig. 5 (a) bond to the 
gold electrodes. Here the Mni2-Ph-Th's easy axis (the z-axis 
in Fig. 5) is approximately perpendicular to the direction of 
the current flow between the electrodes and the Mni2-Ph-Th 
molecule is in its ground state (S= 10 and S z = 10). Spin up 
and spin down are defined relative to the z-axis. 

The calculated spin-resolved transmission probabilities at 
zero source-drain bias are shown in Fig. 12 (a-d) as a func- 
tion of the incident electron energy, E, relative to the Fermi 
energy of the electrodes, There is a very close correspon- 
dence between peaks of the transmission probabilities in Fig. 
12 and the spin-resolved density of states in Fig. 6. The trans- 
mission peaks that are visible at negative energies are due to 
the HOMO and molecular orbitals (MOs) with lower energies 
than the HOMO. Because the MOs close to HOMO are spin 
up electron states, only spin up electrons are transmitted in 
this energy range. However, the LUMO and nearby MOs are 
responsible for peaks at positive energies. As has been dis- 
cussed above, the LUMO is a spin down state, therefore, the 
transmission peak with the lowest positive energy corresponds 
to transmission of only spin down electrons. 

If we regard the gate voltage as an electrostatic potential 
that shifts the energies of MOs rigidly relative to electrodes, 
then by applying negative gate voltages MOs near to HOMO 
are moved closer to electrode Fermi energy. Then at low bias 
voltages only spin up electrons are transmitted, i.e. the current 
through the SMM is purely spin up as is seen in Fig. 12(f). By 
contrast, if a positive gate voltage is applied, then the transmit- 
ted current is spin down at bias voltage below 0.65 eV, Fig. 
12 (e). Therefore our model predicts that Mni2-Ph-Th can be 
used as a spin filter. 

The smooth rise of the current with increasing bias at low 
bias that is seen in Fig. 12 (e) is due to the low energy tail of 
the transmission probability close to LUMO that can be seen 
in the inset of Fig. 12 (b), which is due to the weak hybridiza- 
tion of LUMO, LUMO+1, and LUMO+2 with gold contacts. 
Because the LUMO, LUMO+1, and LUMO+2 have a weak 
presence on both of the ligands that bond to the gold elec- 
trodes (although the amplitude of the LUMO's wavefunction 
on these two ligands is much smaller than that on the out-of- 
plane ligands) (see Fig. 10 (a)) the LUMO, LUMO+1, and 
LUMO+2 are less weakly coupled electronically to the elec- 
trodes than are the other MOs above that LUMO. However, 
since the LUMO and MOs that are close to LUMO in energy 
reside mainly on out-of-plane ligands and the HOMO level 
resides mainly on the magnetic core of the molecule these or- 
bitals are all weakly coupled to the leads and therefore should 
in practice be subject to Coulomb blockade whose main effect 
will be to suppress the tail of the current at low bias in Fig. 12 
(e) and to shift the onsets of the current in Fig. 12 (e) and (f) 
to a somewhat higher values of the bias voltage. As shown in 
the Fig. 10 (c), there are MOs having a strong presence on the 
in-plane ligands, however they are too far above the LUMO in 
energy to affect the transport at the low bias voltages shown 
in Fig. 12 (e) and (f) significantly. 

As we have mentioned above, the strongly spin-polarized 
transmission probabilities that give rise to the spin filtering 
seen in Fig. 12 (e-f) occur because of the strong spin polar- 
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Figure 12. (Color online) (a-d) Calculated spin-resolved transmission probabilities at zero bias and gate voltage as a function of energy E 
relative to Fermi energy £p for a Mni2-Ph-Th molecule bridging a pair of gold electrodes via in-plane ligands in the bonding configuration 
IPBi denned in Fig. 5 (a). Transmission of (a) spin up to spin up, (b) spin down to spin down, (c) spin up to spin down,(d) spin down to 
spin up. Note that the peaks visible in (c) and (d) are a direct consequence of spin-orbit coupling. The insets show the spin up (down) to 
spin up (down) transmission probability at zero bias and gate voltage magnified by a factor of 500 relative to the plots in panels (a ) and (b); 
the tails of the spin down to spin down transmission peaks are responsible for the smooth onset of the spin down current with increasing bias 
in (e). The calculated spin resolved current is also shown for (e) positive and (f) negative gate voltage as a function of bias voltage at zero 
temperature. Solid red ( dashed blue) lines represent spin up (down) current, (e) At low bias voltages for positive gate voltages (gate potential 
at the molecule = +0.2 V) only spin down electrons are transmitted for bias voltages less than 0.65 eV. (f) For the negative gate voltage (gate 
potential at the molecule = —0.2 V) only spin up electrons are transmitted. 



ization of the HOMO and LUMO of the molecule that can be 
seen in Fig. 6. Therefore, qualitatively similar strong spin 
filtering is also predicted for other in-plane bonding geome- 
tries of the molecule and electrodes, provided that the axis of 
spin quantization is always taken to be the easy axis of the 
SMM. However, our model predicts the magnitude of the cur- 
rent passing through the molecule to depend very strongly on 
the bonding geometry. 

In Fig. 13 we present the results of our spin-resolved trans- 
port calculations for a particular bonding geometry in which 
the two out-of-plane ligands labelled OPBi in Fig. 5 (b) bond 
to the gold electrodes. Here the Mn 12 -Ph-Th's easy axis (the 
z-axis in Fig. 5) is approximately parallel to the direction of 
the current flow between the electrodes and the Mni2-Ph-Th 
molecule is in its ground state (S= 10 and S z = 10). Spin up 
and spin down are defined relative to the z-axis. 

There is again a close correspondence between the peaks of 
the transmission probabilities in Fig. 13 (a-d) and the spin- 
resolved density of states in Fig. 6. 

The transmission peaks at low negative energy, due to the 
HOMO and molecular orbitals with lower energies than the 
HOMO, for IPBx and OPBi are similar (see Fig. 12 (a-d) 
and Fig. 13 (a-d)). That is, again only spin up electrons are 
transmitted in this energy range (see also Fig. 13 (f)). 



However, in the OPBi case in Fig. 13 (a) there is a long 
low energy tail of the spin up to spin up transmission proba- 
bility that extends far below the LUMO in energy. This tail is 
due to the strong level broadening that results from the strong 
hybridization of the LUMO+5, LUMO+14, and LUMO+17 
with the gold contacts. These MOs are indicated by the solid 
circles in the projected density of states (PDOS) plot in Fig. 
10. The hybridization is strong because these orbitals reside 
mainly on the out-of-plane ligands that are bonded to the elec- 
trodes. One of these orbitals is displayed in Fig. 10 (b). Even 
though the LUMO is a spin down state, the off -resonant tun- 
neling through the molecule mediated by this hybridization is 
responsible for the spin up current at low bias voltage in Fig. 
13 (e) since the LUMO+5, LUMO+14, and LUMO+17 are 
spin up orbitals as is seen in Fig. 10. Note that the LUMO+5, 
LUMO+14, and LUMO+17 have much more weight on the 
ligands that connect to the gold in the OPBi bonding geom- 
etry than does the LUMO. For this reason the hybridization 
of the spin up LUMO+5, LUMO+14, and LUMO+17 with 
the gold is much stronger than that of the spin down LUMO 
and hence the broadening of the LUMO+5, LUMO+14, and 
LUMO+17 due to their coupling with the gold is also much 
stronger. Hence the spin up current at low bias voltages in 
Fig. 13 (e) is much stronger than the spin down current. 
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Figure 13. (Color online) (a-d) Calculated spin-resolved transmission probabilities at zero bias and gate voltage as a function of energy E 
relative to Fermi energy Ep for a Mni2-Ph-Th molecule bridging a pair of gold electrodes via out-of-plane ligands in the bonding configuration 
OPBi defined in Fig. 5 (a). Transmission of (a) spin up to spin up, (b) spin down to spin down, (c) spin up to spin down,(d) spin down to spin 
up. Note that the peaks visible in (c) and (d) are a direct consequence of spin-orbit coupling. Calculated spin resolved current for (e) positive 
and (f) negative gate voltage as a function of bias voltage at zero temperature. Solid red (dashed blue) lines represent spin up (down) current, 
(e) For positive gate voltages (gate potential at the molecule = +0.2 V) at low bias voltage (0 <Yuas <0.53 V) most of the transmitted 
electrons are spin up. For Vuas >0.53 V the current is not spin polarized, (f) For negative gate voltages (gate potential at the molecule = —0.2 
V) at low bias voltage only spin up electrons are transmitted. 



The very strong hybridization of the LUMO+5, LUMO+14, 
and LUMO+17 with the gold electrodes in the the OPBi 
case and the resulting level broadening is partly responsi- 
ble for the low bias tail of the current in Fig. 13 (e) being 
roughly two orders of magnitude stronger than that in Fig. 
12 (e). For the same reason, we predict tunneling associated 
with the LUMO+5, LUMO+14, and LUMO+17 levels for the 
molecule in the bonding geometry that we are considering in 
Fig. 13 to be immune to suppression by charging effects, i.e., 
there should be no Coulomb blockade of transport associated 
with the LUMO+5, LUMO+14, and LUMO+17. By contrast, 
since the HOMO and nearby levels reside mainly on the mag- 
netic core of the molecule they are much more weakly coupled 
to the leads and much less broadened. This is evident from the 
lack of visible tails of the resonant peaks at negative energies 
in Fig. 13 (a) and the correspondingly very abrupt character 
of the current steps in Fig. 13 (f). Therefore transport via 
the HOMO and nearby levels should be subject to Coulomb 
blockade whose main effect will be to shift the onset of the 
current in Fig. 13 (f) to a somewhat higher value of the bias 
voltage. 

If we regard the gate voltage as an electrostatic potential 
that shifts the energies of MOs rigidly relative to electrodes, 
then applying negative gate voltages shifts MOs that are near 
to HOMO closer to the electrode Fermi energy. Then at low 
bias voltages purely spin up electrons are transmitted, simi- 



larly to the case of IPBi bonding to the electrodes. I.e. the 
current through the SMM is purely spin up as is seen in Fig. 
13 (f). But, in contrast to IPBi, if a positive gate voltage is 
applied, then the transmitted current at low bias voltage (0 
< Vbias < 0.53 V) is predominantly spin up, Fig. 13 (e). 
However, at bias voltages above 0.53 V, when the LUMO 
and nearby MOs which have spin down become available for 
transport, the transmitted current ceases to be spin polarized. 

Therefore our model predicts that Mni2-Ph-Th can be used 
as a spin filter. However, we predict that it should exhibit 
differing behaviour depending on the orientation of molecule 
relative to the electrodes. At negative gate voltages for both 
in-plane bonding and out-of-plane bonding configurations we 
predict that the molecule should transmit spin up electrons. 
However, for positive gave voltages at low bias voltage, in 
the case of in-plane bonding geometries we predict spin down 
electrons to be transmitted, whereas for out-of-plane bonding 
we predict spin up electrons to be transmitted. For positive 
gave voltages at higher bias voltages the transmitted current 
ceases to be spin polarized. 

Note that although the transport calculations presented in 
this work do not include the effects of Coulomb blockade ex- 
plicitly whereas transport via the molecular HOMO and near 
HOMO levels is expected to be subject to Coulomb blockade, 
our prediction of effective spin filtering occurring in trans- 
port via the molecular HOMO and near HOMO levels is ex- 



14 



pected to be valid in the presence of Coulomb blockade at 
the Coulomb blockade threshold. This is because of at the 
Coulomb blockade threshold the only available state for elec- 
tron to pass through the molecule is the HOMO or a near 
HOMO orbital all of which are predicted to have spin aligned 
with easy axis 25 (see Fig. 4 (a,c)). Therefore spin flip pro- 
cesses (that would degrade the spin filtering) are not permitted 
energetically on the molecule during transport via the HOMO 
and near HOMO levels even for incoherent transport in the 
Coulomb blockade regime. 



V. CONCLUSIONS 

In this paper we have introduced a tight-binding model of 
single molecule magnets that incorporates both the spin and 
spatial aspects of the molecular electronic structure. It is able 
to provide a realistic description of the known properties of 
Mni2 single molecule magnets for both the neutral and neg- 
atively charged molecules. The properties calculated within 
the present model that are consistent with the available ex- 
perimental data include the total spin of the molecule, the 
spins of the individual Mn ions, the magnitude of the mag- 
netic anisotropy barrier and its orientation and the size of the 
HOMO-LUMO gap. The spins and spatial locations in the 
molecule of the HOMOs predicted by the model are consis- 
tent with the results of density functional theory based cal- 
culations available in the literature. 2 1_23 ' 26 ' 28 Different density 
functional theory based calculations have yielded differing re- 
sults for the spins and spatial locations of the LUMOs. 22 ' 23 
The predictions of the present model for the spins and spatial 
locations of the LUMOs are consistent with experimental data 
for the corresponding negatively charged molecules. 

The description of the properties of the Mni2 family of sin- 
gle molecule magnets provided by the model is comparable to 
and in some respects more realistic than that provided by den- 
sity functional theory that has been improved by the inclusion 
of a phenomenological Hubbard U parameter. However, cal- 
culations based on the present model are much simpler and 
require far less computational resources and compute time 
than density functional theory-based calculations. We have 
therefore been able to investigate the electronic and spintronic 
structure and transport properties of more complete models of 
Mni2 single molecule magnets than have been studied previ- 
ously. In particular, we have presented the results of the first 
transport calculations for Mni2 single molecule magnets with 
complete sets of ligands none of which have been truncated or 
replaced with hydrogen atoms. 

In the models of transport in Mni2 single molecule mag- 
nets that have been studied previously theoretically, the lig- 
ands have always behaved as simple tunnel barriers between 
the magnetic core of the molecule and the electrodes and thus 
transport has been in the Coulomb blockade regime. Here we 
have shown that for Mni2 single molecule magnets covered 
with thiolated or methylthio terminated aromatic benzoate lig- 
ands, the situation can be very different: We predict that for 
these molecules the molecular LUMO and molecular orbitals 
close in energy to the LUMO reside on some of the ligands. 



This has important novel implications for the transport prop- 
erties of these single-molecule magnets: In the case where a 
molecular orbital resides on one or more ligands that bond to 
an electrode, we predict low temperature conduction via that 
molecular orbital to be in the coherent resonant or off-reso nant 
quantum regime and not subject to Coulomb blockade that has 
until now been assumed to always govern electrical conduc- 
tion through Mni2 single-molecule magnets. Because we find 
the LUMO and nearby orbitals to be localized on ligands that 
are aligned approximately parallel to the magnetic easy axis 
of the Mni2 single molecule magnet, we propose that observa- 
tion of the absence of Coulomb blockade at low source-drain 
bias and positive gate voltages can be used as an experimental 
signature identifying specific realizations of single molecule 
magnet transistors in which the magnetic easy axis and mag- 
netic dipole moment are aligned approximately parallel to the 
direction of current flow through the molecule. We also pre- 
dict that Mni2 single molecule magnets covered with thiolated 
aromatic benzoate ligands should behave as effective spin fil- 
ters, and that this should be the case even for transport at low 
bias via the molecular orbitals that reside primarily on the lig- 
ands and not on the magnetic core of the molecule. While 
we have predicted that the LUMO should lie on the ligands of 
Mni2-benzoate and its derivatives, it is reasonable to expect 
the LUMO and/or the HOMO to lie on the ligands of some 
other Mni2 single molecule magnets as well. Examples of 
such systems may be Mni2 with ligands derived from poly- 
acetylene and polythiophene that have small HOMO-LUMO 
gaps, 1.4 eV 76 and 0.85 eV 77 , respectively. 

If these predictions are confirmed experimentally, the im- 
plications for the field of transport in single molecule magnets 
will be significant: Transport through Mni2 single molecule 
magnets that is not in the Coulomb blockade regime will have 
become experimentally accessible. Also, single molecule 
magnet transistors with a known orientation of the magnetic 
easy axis (and hence a known orientation of the spins carried 
by the spin-filtered current) relative to the direction of current 
flow through the single molecule magnet will have been real- 
ized. 
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Appendix A 

In this Appendix we evaluate the matrix elements of the 
spin-orbit coupling Hamiltonian H so between valence or- 
bitals i and i' of atoms (3 and /?' with spin s and ,s'. Eq. 4 
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of Section II C in matrix form is 

1 1 dV a (\r-r a \) 

?.m z c z r — r a d\r — r Q 

(Al) 

where 'fj^ is the i th atomic orbital of the /3 th atom with spin 
s. The largest terms in the summation in Eq. Al are those 
in which (atom) a is the same as (5 or j3' or both. In what 



follows, we shall consider only these terms. 

Then for f3 = /?', the spin-orbit coupling Hamiltonian ma- 
trix element is due to intra-atomic spin-orbit coupling. We 
will refer to it as E^ 3 ,. a = Ef sB .y s , B . The atomic orbital 
wave function ^i S p can be expressed as the product of a radial 
wave function Rr,u an d directed atomic orbital \f3,l i ,d i ,s). 
Here k is the angular momentum quantum number, di may be 
s,Px,Py,Pz, d xy ,d xz , ... depending on the value of k, and s 
is the spin quantum number. Then 



= (P,li,di,8\8.Lp\M,di.y)(RuA ^ |r ^ (A2) 

The S ■ operator does not change the angular momentum quantum number ly, therefore E^^ s ,.g is non-zero only when 
li = lit. Explicit expressions for the resulting matrix elements of S • in Eq. A2 have been given in Ref. 58. The radial 

integrals, (Rp,i t 2 m' A c 1 | r-r 3 | rf <j(|r-~rg'|j" > \Rp<k) = e /3 'i' are me s P m -orbit coupling constants which are available in the literature 
for various atoms and ions. The numerical values of these constants are discussed in Sec. III. 

For the z'«fe r-atomic contributions to the spin-orbit coupling, in Eq. Al j3 ^ 0' and the summation over a reduces to two 
dominant terms, those for which a = /3' or a = (3. The term a = j3' is 

^^m^lr-r^l d\r-r ,\ 8 1 VI*W/»'> (A3) 

= Yt^W^W > <PW I o~~~2~lz i 1 dV ^ r ~ r ^ S • Lp ^p,) (A4) 

2m 2 c 2 |r — r g >\ d\r — v rj i\ 



The summation in expression (A4) is over all of the atomic 
orbitals of the f3'th atom that form a complete orthogonal 
set. However, since our goal is to construct a generalization 
of extended Hiickel theory that includes spin-orbit coupling 
we retain in the sum only the atomic valence orbitals. Then 
the expression (A4) reduces to J2j ^iP;j0'-^js?s'-0' where 
Dip.jp, = (%^^}S s , sr 

The remaining term a = (3 of Eq. (Al) is 



_J 1 dVp(\r-r p \) 

2m 2 c 2 |r — rg d\r — rp\ 



S-L/jIWiv^/) 



i's'ft' | 



1 dVp(\T-T P \) 



S-Lfll*^)* 



2m 2 c 2 |r — rp\ d\r — ra| 
since S • hp is Hermitian. The last expression is evaluated in a 



similar way to Eq. (A3). Finally, collecting the above results, 
Eq. (Al) for the matrix elements of the spin-orbit coupling 
Hamiltonian reduces to 



< iS /3|£HzV/3'} ~ E^ s ,, 8 Sp 0l + (1 - <V) 



ZT'intxa 

iPi3p ,£j jai'8';P' (A5) 



+ [Di'p'-j/3Epl s . B }*) 

Here the first term on the right-hand side is the intra-atomic 
contribution and the remaining terms are the inter-atomic con- 
tribution. 
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